Nonequilibrium electron spectroscopy of Luttinger liquids 
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Understanding the effects of nonequilibrium on strongly interacting quantum systems is a challenging prob- 
lem in condensed matter physics. In dimensions greater than one, interacting electrons can often be understood 
within Fermi-liquid theory where low-energy excitations are weakly interacting quasiparticles. On the contrary, 
electrons in one dimension are known to form a strongly-correlated phase of matter called a Luttinger liquid 
(LL), whose low-energy excitations are collective density waves, or plasmons, of the electron gas. Here we 
show that spectroscopy of locally injected high-energy electrons can be used to probe energy relaxation in the 
presence of such strong correlations. For detection energies near the injection energy, the electron distribution is 
described by a power law whose exponent depends in a continuous way on the Luttinger parameter, and energy 
relaxation can be attributed to plasmon emission. For a chiral LL as realized at the edge of a fractional quan- 
tum Hall state, the distribution function grows linearly with the distance to the injection energy, independent of 
filling fraction. 



Over the last decade, experimental advances in nanostruc- 
ture fabrication have brought a resurgence of interest in the 
LL model because of the possibility to test its peculiar pre- 
dictions Qiiii. 

Defining signatures of a LL such as 
spin-charge separation 001, charge fractionalization 11,1^, 
and the power-law suppression of the local electron tunnel- 
ing density of states fioi H, 12, 13, 3 have been experi- 
mentally verified. Recently, LLs driven far from equilibrium 
have begun to receive attention |15 , E In, 18]. Studying 
these systems offers the possibility to characterize novel as- 
pects of electron-electron interactions and to understand en- 
ergy relaxation processes that have not been apparent in the 
above-mentioned equilibrium experiments. 

Here we consider a LL driven out of equilibrium by local 
injection of high-energy electrons, far away from any con- 
tacts, at a fixed energy. Their spectral properties are extracted 
at another spatial point some distance away by evaluating the 
average tunneling current from the LL into an empty reso- 
nant level with tunable energy. In this work, we consider both 
standard (non-chiral) and chiral LLs, which are realized at the 
edge of fractional quantum Hall systems |3, 4, 12, 13, 14]- 

For the standard LL and for probe energies slightly below 
the injection energy, we find that the inelastic component of 
the current shows a power law behavior as a function of the 
difference between injection and detection energy, with an ex- 
ponent that continuously evolves as the interaction parameter 
is varied. We develop a perturbative approach which shows 
how injected electrons can relax by emitting plasmons inside 
the wire. 

For a chiral LL at the edge of a fractional quantum Hall state 
from the Laughlin sequence, an essentially exact calculation 
of the tunneling current is possible. Here, the inelastic part of 
the electron cuiTent increases in a linear fashion as the probe 
energy is lowered from the injection energy towards the chem- 
ical potential of the edge state, despite a decreasing tunneling 
density of states for electrons. This behavior is compatible 
with our result for the standard LL in the limit of strongly re- 
pulsive interactions. For probe energies close to the chemical 
potential, the chiral LL is far from equilibrium: the electron 



x=0 Luttinger liquid x=L 

FIG. 1: The proposed experimental setup. Hot electrons are in- 
jected from the source resonant level at x = 0, and are collected at 
the probe resonant level at x = L. System parameters are set (see 
text) so that the source (probe) occupancy is fixed to be full (empty). 
Spectral properties of the injected electrons are extracted by measur- 
ing the tunneling current between the edge and the probe (indicated 
by the arrow). 



spectral function approaches a finite value, in striking contrast 
to the power law decrease towards zero in equilibrium. In ad- 
dition to the inelastic contribution to the probe current, in a 
chiral LL there always is an elastic contribution, indicating 
that a finite fraction of electrons travels from the injection to 
the probe site without loosing energy. 

Electrons with charge eo are injected into the LL from a 
resonant level (source) with energy Ei = cqVi > at position 
X - Q (see FigQJ. Energy relaxation is studied by coupling a 
second resonant level (probe) with energy E2 = eoVa > to 
the LL at position x - L (downstream for the chiral LL), and 
by computing the tunneling electron current between the LL 
and that level. The two levels are coupled to the LL via tun- 
neling amplitudes 771 and 772, respectively. In addition, source 
and probe dots are coupled to reservoirs held at chemical po- 
tentials yUi and yU2 via tunneling amplitudes Ai and Ai- The 
chemical potential of the LL is taken to be zero. We assume 
the level broadening due to tunnel couplings to be small in 
comparison to both E\ and Ei, and therefore consider the cur- 
rent in the sequential-tunneling regime. Further, we assume 
Ai » 771 with //I > £■] so that the source occupancy is con- 
strained to one, and A2 » 772 with 112 < £2 so that the probe 
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FIG. 2: TXinneling current for the standard case for various in- 
teraction. The current is plotted for zero temperature and includes 
the leading contribution in AE/Ei. The inelastic contribution for 
E2 < El shows a power law decay as a function of increasing AE 
with an exponent that depends on the interaction parameter (see text). 
A level broadening of 0.01£i is used for the elastic peak. 



FIG. 3: Tunneling current for the chiral case at various tempera- 
tures. The current shows an elastic contribution at £1 = E2, and an 
inelastic contribution for E2 < Ei which increases as energy trans- 
fer is increased. The broadening of the elastic peak is included as in 
the standard case. The same level broadening as the standard case is 
used here. 



occupancy is fixed at zero. 

We first focus on the standard LL and consider spinless 
electrons, for which the interaction strength is described by 
a single parameter K |2]. The case K - 1 describes non- 
interacting electrons, K < I corresponds to repulsive inter- 
actions, and K > I to attractive ones. We use the non- 
equilibrium Keldysh formalism lfl9ll to calculate the current 
flowing into the probe dot to leading order in the tunneling 
amplitudes 771 and 772, details are described in the appendix. 
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Here, r(x) is the gamma-function, u is the velocity of plasmon 
excitations, a denotes the short distance cutoff of the theory, 
and y - K(l/K - l)'^/4 > 0. When calculating the current 
Eq. ([T]i, the limit of large interdot separation u/L <K E\,E2 
was taken. In the non-interacting limit {y 0) the quantity 
in the square brackets is a representation of the delta-function, 
and Eq. ^ reduces to / oc 5{AE). When the interactions are 
turned on, the elastic peak gradually broadens to give rise to 
an inelastic contribution which shows a power law decay as 
a function of increasing AE, with an exponent that continu- 
ously evolves as a function of the interaction parameter. For 
strong enough interactions with y > 1/2, the elastic peak van- 
ishes and the remaining inelastic contribution monotonically 
increases with a power law which again evolves as a function 
of the interaction parameter. The result Eq. ([T) is plotted in 
Fig|2] Broadening of the peak is included in the figure to re- 
flect the finite width of the resonant levels due to the couplings 
to the reservoirs and the wire. 

In the limit of weak interactions with Luttinger parameter 
K close to one, energy relaxation as described by Eq. ^ can 
be interpreted by using lowest order perturbation theory in the 
interaction strength. Interactions can be decomposed into for- 



ward scattering between electrons near the same Fermi point 
with amplitude gn, and between electrons near opposite Fermi 
points with amplitude g2- The g4-interaction merely renor- 
malizes the fermion and plasmon velocities and cannot give 
rise to relaxation. In a spatially homogeneous LL, the g2- 
process cannot give rise to energy relaxation either due to the 
simultaneous requirement of momentum and energy conser- 
vation. However, because of the local nature of injection and 
collection processes considered here, an electron is capable 
of exploring virtual momentum states in connection with tun- 
neling, and a consecutive inelastic process can both conserve 
momentum and produce a final state with the same total en- 
ergy as the source state. Here, we consider the lowest order 
inelastic process proportional to y oc ^ero temperature, in 
which an electron is transported from the source to the probe 
while emitting a single plasmon inside the wire. First, the 
electron in the source tunnels into a right moving momentum 
eigenstate whose energy may be different from E\ . In the sec- 
ond step, a left moving plasmon with energy AE is emitted 
via a g2-process, and the right moving electron is scattered 
into another wire state such that momentum is conserved. Af- 
ter propagation along the wire, the electron tunnels into the 
probe. Alternatively, tunneling into the wire can be elastic, 
and the plasmon can be emitted when tunneling into the probe. 
One finds that the matrix element for these processes scales as 
1 / yl\AE\. This can be understood by multiplying the matrix 
elementfor plasmon emission, which increases as yl\AE\, with 
the time available for plasmon emission, which diminishes as 
1/IAiil due to the energy-time uncertainty principle. The tun- 
nel current can then be computed using Fermi's golden rule, 
and correctly reproduces the inelastic component / oc y/A£ 
of equation ([TJ to order y. 

Next, we consider tunneling into a chiral LL at the edge of a 
fractional quantum Hall state from the Laughlin sequence. We 
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focus on the filling fraction v - 1/3, where the area occupied 
by one electron is threaded by three quanta of magnetic flux. 
The calculation of the steady state current proceeds along the 



same lines as for Eq. ([T]), with the difference that we were able 
to obtain an exact expression for all values of AE and for finite 
temperature, 
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Here, AX - X\ - X2, X, = Ei/kBT. At zero temperature, 
the expression for the current simplifies to / oc E'^S{AE) + 
6e{AE){E'\ + El)AE where AE ^ Ei - E2. The cuiTent, plot- 
ted for both zero and finite temperatures in Fig|3] has two main 
contributions; elastic and inelastic. The peak is due to elec- 
trons that were elastically transported from the source to the 
probe. Second, there is a broad inelastic contribution that ex- 
tends over the range E2 < Ex, and that grows monotonically 
as E2 is lowered. For E2 < E\, the current increases linearly 
with AE. We have confirmed that a similar inelastic contri- 
bution to the current is also present for the Laughlin filling 
fraction v - 1/5. In this case, an exact computation at zero 
temperature shows again that /;„£! AE for E2 < E\. This 
suggests that the linear upturn in the current below E\ may be 
a generic feature at all Laughlin filling fractions. 

For a non-interacting chiral Fermi liquid, which describes 
the edge excitations of an integer quantum Hall state, hot elec- 
trons do not relax. In addition, the weight of the elastic peak is 
reduced as the temperature is increased. This reduction is due 
to Pauli blocking of states by thermally excited edge electrons 
residing above the chemical potential. When interactions are 
present, Fig[3] shows an overall increase in the elastic current 
with temperature. This reflects the increase in the tunneling 
density of states with temperature and constitutes a clear sig- 
nature of LL physics. 

The setup of Fig [1] is ideal for directly extracting the elec- 
tron energy distribution, /(£), and spectral function, A(E), 
inside the wire at a spatial point far from the injection site. 
With the probe occupancy constrained to be empty, the tun- 
neling current is given by /empty = 'eo|772pG*(£'), while 
a similar evaluation with probe occupation held full gives 
/full - ieo\ri2?'G^(E) [20]. Once the two currents are ob- 
tained, both f{E) and A{E) can be extracted by expressing 
the lesser and greater Green functions, C^iE) - if(E)A{E) 
and G^{E) - - f(E))A{E), in terms of electron distribu- 
tion function and spectral weight. At zero temperature and for 
V = 1/3, f(E2) and A(E2) vahd for < £2 < £1 read 
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where E* = yj6n\rii\^a^E^Ju^h^ separates two energy 
regimes. We note that E* can be parametrically larger than 
the level widths such that our sequential tunneling approxima- 
tion stays valid. In the high-energy regime and for small en- 
ergy transfers (E* <s £2 < £1), f(E2) ~ l2n\rii\^a^AEIw'rr', 
which shows that the linear upturn in the current below £1 is 
also reflected in the distribution function. In the same regime, 
we find that the spectral function does not deviate strongly 
from its equilibrium expression (with rji - 0). In the low- 
energy regime (0 < £2 < £*), f{E2) smoothly approaches 
one and the spectral function approaches a finite value. The 
latter is in stark contrast to the equilibrium case. 
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Appendix: We now provide the theoretical basis for the 
derivation of Eqs. ( 1112b . The system is modeled by the Hamil- 
tonian H - Hu^ + //dot + //tun, where //ll models the LL, 
//dot - E\il/\il/\ + E2^2'^2 the two resonant states, and //tun 
describes the tunneling of electrons between the wire and 
the two resonant levels, i/'i {^2) are electron operators of 
the source (probe) with occupation numbers (iAJi/'i) - 1 and 
{'^['^2) - 0. The standard LL Hamiltonian reads 101 



//ll=— r 



(5) 



where K is the LL parameter, and the left and right moving 
boson operators satisfy [(p]i(x),(pR(x')] - -[(Pl(x),(Pl(x')] = 
inK sgn(x - x'). One-dimensional electron densities are given 
by Pr,l{x) - {d_x(fiR,Lix))l1jT: and u denotes the plasmon veloc- 
ity. To simplify the notation, we use the units where Ti - \ and 
ks - I- The tunneling Hamiltonian is given by 



//tun = T]ilJ/ilp'(x = 0) H- T]2l//2lp'(x = L) -H h.C. 



(6) 



where iff(x) - i//r(x) + if/iix)- The electron operators can 
be bosonized as iffRx(x) - exp[i(K±cf>R(x) + K^<pi(x))]/ y/lna 
with K± = (K^^ + l)/2. The expectation value of the current 
reads 



/ = <rc{/ci(fi)e-'X''"^-«))o, 



(7) 
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where all operators are written in the interaction picture 
with respect to //ll + ^^dot- The current is computed using 
the nonequilibrium Keldysh technique |19], and indicates 
time-ordering of the operators on the time-loop contour c. The 
"classical" component of the current operator is the symmetric 
combination of the operator on the forward {+) and backward 



(-) parts of the Keldysh contour, i.e. /d(0 - (h{t) + I-(t))/2, 
where I±(t) — -;eo[//-t, i/'j ^i^2,±]- Upon imposing the con- 
straints on the resonant level occupancies and taking the limit 
of large inter-dot separation, the time of propagation L/u 
drops out and we arrive at the following expression for the 
steady state current to leading order in 771 and 772, 



J 
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Here, 7 = K^K, y' = (KI+KJ)K, and t-.j = f,-fj. The ordering 
on the Keldysh contour for a related problem is described in 
112 in . The factors of correlation functions. 
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can be interpreted as the tunneling in and out density of states, 
and the H-matrices, 



^ G^(-f3)G-(f24)' 



(10) 



describe the propagation of electrons along the wire. In prin- 
ciple, the current in equation (|7]) contains another term pro- 
portional to only \r]2\^ that describes the tunnehng of electrons 
into the probe from thermal excitations in the wire. However, 
for low temperatures (T <K £1, £2), this contribution is expo- 
nentially suppressed. 

The current for the chiral LL Eq. (|2| can be derived in a 
similar fashion. First, we note that the Hamiltonian is analo- 
gous to that in Eq. (|5]l, but with only one boson field, say <pR, 
and with K replaced by the filling fraction v. The tunneling 
Hamiltonian is identical to equation (|6]l, and the electron op- 
erator is now bosonized as i//(x) - yflna. The formal 
expression for the current is still given by equation (|7|l- Using 
similar steps as above, one finds 

dt2dhdhe-''^-'-^''^'''' 

CO 

X {iG<i^{-t2)){iG\i^ihM^t"y - ^tp- (11) 

The correlation functions and the D-matrices are again given 
by equations ( 19110b . 
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